Space symmetries draw elasticity theory 
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The foundation of continuum elasticity theory is based on two general principles: (i) the force felt 
by a small volume element from its surrounding acts only through its surface (the Cauchy principle, 
justified by the fact that interactions are of short range and are therefore localized at the boundary); 
(ii) the stress tensor must be symmetric in order to prevent spontaneous rotation of the material 
points. These two requirements are presented to be necessary in classical textbooks on elasticity 
theory. By using only basic spatial symmetries it is shown that elastodynamics equations can be 
derived, for high symmetry crystals (the typical case considered in most textbooks), without evoking 
any of the two above physical principles. 
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I. INTRODUCTION 

When teaching elasticity theory one is faced with the 
question of whether it is necessary to introduce series of 
postulates (to be recalled below) in order to arrive to the 
final elastodynamics equations (like Hookes law. Lame 
equation, etc.), or is there a simpler way bypassing gen- 
eral principles (except using the Newton law). Moreover, 
traditionally, very quickly the notion of tensors (of second 
and fourth rank, see below) is used, which adds new ob- 
stacles to many students. Here we give a direct derivation 
of elasticity equations by evoking only elementary space 
symmetries. The notion of stress tensor (of second rank) 
appears only at the very end as a natural consequence 
of elementary symmetries, and not as an initial notion 
entering the ingredients of the theory. The main spirit of 
this paper is presented in section IL 

Let us first recall the traditional foundation principles 
of elasticity theory. In classical textbooks, such as those 
of Landau and Lifchitz^, Fcynman^, and others^"—, this 
theory is based on the following three steps. 

(i) The force felt by a given material element acts only 
on the surface. Consider a small volume element of the 
material having a volume AT^, enclosed by a surface AA. 
After deformation that element would feel a force distri- 
bution from the surrounding medium acting only through 
its surface. The elementary force di acting on a portion 
dA of the closed surface AA has the form 



di = andA, 



(1) 



where cr is the stress tensor and n is the normal unit 
vector. That the force is only localized at the surface 
(and not to the bulk element) is often justified by the fact 
that the interactions are of short range. For example, it is 
stated in the book of Landau and Lifchitz^ (p. 4): We can 
say that the forces which cause the internal stresses are, 
as regards the theory of elasticity, "near-action" forces 
[. . .] Hence it follows that the forces exerted on any part 
of the body by surrounding parts act only on the surface 
of this part. 



The fact that the force must be localized at the sur- 
face and have the form ([T|) is known as the Cauchy 
principle. The resultant force applied on the volume 
element is obtained upon integration over the surface. 



Iaa '^^d,^ 



J^y V • crdV, where we have used the 
divergence theorem. This means that the force per unit 
volume is simply given by V • cr. From this and the New- 
ton law we can write 



pu 



(2) 



where p is the mass density, u is the displacement vec- 
tor, and pii represents the inertial force per unit volume 
(double dots designate the second time derivative). 

In the book of Feynman? the fact that the stress ten- 
sor acts on the surface is justified as follows (p. 31-6 of 
Ref. |2|) : Consider a body of some elastic material -say a 
block of jello. If we make a cut through a block, the mate- 
rial on each side of the cut will, in general, get displaced 
by internal forces. Before the cut was made , there must 
have been forces between the two parts of the block that 
kept the material in place; we can define the stresses in 
term of these forces. 

(ii) The stress tensor is symmetric. The stress tensor 
a must be symmetric, i.e. 

aij = aji. (3) 

The justification of this condition varies from one book 
to the other. In Ref. [l|, calculating the torque of the 
elastic forces on the volume element, it is stated that if 
the stress tensor is not symmetric then the total torque 
would have not only a surface contribution, but also a 
bulk contribution, which would contradict the fact that 
the forces are localized to the surface only. The bulk 
contribution vanishes if the stress tensor is assumed to 
be symmetric. A variant of this proof can be found in 
the book of Feynman^ where it is said (p. 31-11): We 
can also show that Sij is a symmetric [stress] tensor by 
looking at the forces on a little cube of material. [. . .] Now 
there must be no torque on the cube, or it would start 



spinning. Similar arguments are given in the books of 
Love^, Timoshenko and Goodier^, and Kittel'^. 

(iii) The stress tensor has a linear relationship with 
the strain tensor. The stress tensor has a proportionahty 
relation with the strain tensor e^j = ^{diUj+djUi), where 
Ui is the i-th cartesian component^ of the displacement 
vector u. For an isotropic material, we have 



aij = XekkSij + 2 fie 



-i], 



(4) 



where A and fi are the two Lame coefficients.^ In most 
books (see for example Ref. |2|, p. 31-7) this is derived 
from the energy by saying that the elastic energy is a 
quadratic function of the strain (F — h^ijkiUj^ki), by 
(i) and (ii) evoking the symmetry of the stress tensor, 
and by using some general mathematical relations obeyed 
by the rank four tensor Kijki for an isotropic medium. 
Remark in passing, that our own experience is that many 
students are quite reluctant already regarding stress and 
strain tensors (second-rank tensors), and this becomes 
even worse when a fourth-rank tensor is evoked. 

Reporting @ into ^ one arrives to the elastodynam- 
ics equations for an isotropic material. 



pu 



/iV^u 



(A + /i)V(V-u). 



(5) 



The objective of this article is to show that for high 
symmetry crystals (such as an isotropic and cubic crys- 
tals in 3D and isotropic, square and hexagonal in 2D 
crystals) none of the above three points is necessary in 
order to arrive to the elastodynamics equations. Rather, 
evoking only elementary space symmetry is sufficient. 



II. DERIVATION OF ELASTICITY EQUATIONS 
FROM ELEMENTARY SPACE SYMMETRIES 

In order to make the derivation as pedagogical as possi- 
ble, we shall consider first a two dimensional system. Let 
u(a;, y) be the displacement vector having (u, v) as Carte- 
sian components along the x and y directions. An elastic 
force takes place only if \i{x,y) is non constant, other- 
wise all the material element would displace by the same 
quantity, corresponding to a rigid translation. The fact 
that u(x, y) is non constant is expressed by non zero first 
derivatives, and possibly higher order derivatives. For the 
sake of abbreviation, we shall denote the first derivatives 
as Ux = du/dx, Uy = du/dy, and so on. Second deriva- 
tives will be denoted as u^y = d'^u/dxdy, v^x = d'^vjdx^ , 
and so on. Thus the forces are functions of derivatives of 
u and V. 

Like in most physical theories (elecromagnetism, dif- 
fusion, hydrodynamics,®-...) we shall seek for the lowest 
derivatives that provide relevant results. It will turn out 
that we need to expand the force up to second deriva- 
tives of u(x,y). The force (with Cartesian components 
denoted as U and V) acting on a material element is a 
function of derivatives of u(a;, y), U = U{ui,Vi,Uij,Vij) 
and V 



could also be considered if one is interested in the so- 
called "second gradient" theory, a step beyond the clas- 
sical elasticity theory, a question which is beyond our 
objective. In two dimensions there are four distinct first 
derivatives Ux,Uy,Vx, Vy and six distinct second deriva- 
tives Uxx,Uxy,Uyy,Vxx,Vyy,Vxy From Ncwtou law, the 
acceleration is proportional to the force (we consider a 
force per unit volume for the sake of comparison with 
traditional theory) , and the two components of displace- 
ment field obey 



pil — U, pv = V. 



(6) 



Since we are interested in linear elasticity, the force 
components U and V are developed in a Taylor expansion 
in terms of their arguments by keeping only linear terms 
in u and v. Considering the terms with first derivatives, 
we have explicitly 



pu 
pv 



du 



b'l 



y 
d'v,, 



(7) 



where the coefficients a,h,... are real constant. Let us 
start with the case of an isotropic material. The evolution 
equations must be separately invariant under each of the 
following symmetries: 



X — ?> —X and u — >■ —u, 
y — > — y and v — >■ —v. 



(8) 



This operation corresponds to mirror symmetries, and 
the governing equations should be invariant under this 
operation. This invariance is valid for any crystal that 
has mirror symmetries along x and y directions. This 
is in particular the case for square and hexagonal (six- 
fold) symmetries. For example, using the first of this 
symmetry transforms Eqs. ([7]) into 



pu 



bu,, 



dv,. 



pv — +a Ux ^ b Uy — c Vx + d Vy. 



(9) 



The equations must remain invariant, which implies that 
a = d = b' = c' = 0. Using the second symmetry ^ 
leads to the conclusion that b = c = a' = d' = 0. 

In conclusion, no term with first derivatives must be 
present in the equations ([7]). We have thus to consider 
now second derivatives. The Taylor expansion of the 
functions U and V yields 

pil = aUxx + PUyy + "fVxy + aUxy + /3Vxx + IVyy 

pi = a'Vxx + P'Vyy + j'Uxy + Oi V xy ^ [3' Uxx + l' Uyy 

(10) 

where the coefficients a, /?, ... are real constants. As we 
did with first derivatives above, using the symmetry op- 
erations ((HI) leads quite simply to the conclusion that all 
coefficients having tilde must be zero, getting 



V{u 



2 , 1^2 , ^ij , iJxj 



Vij). Higher order derivatives 



pu 
pv - 



■ O^Uxx 

a'vxx 



PUyy +JVxy, 
- (3'Vyy +j'Uxy 



(11) 



As stated above, the symmetry operation ([8]) is valid 
for any crystal symmetry that enjoys the mirror sym- 
metries ([8]). This includes, in two dimensions, isotropic 
materials, square symmetry and hexagonal crystals. 

We have not yet exploited all symmetry operations. 
For example, for an isotropic material equations (jlip 
must be invariant upon rotation of the coordinate system 
with an arbitrary angle 9 around the 2-axis. Let x' and 
y' designate the new coordinates and u' and v' the new 
displacement components after rotation, we then have 



a;' = xcos6' — y sin6', y' 



X sm ( 



y cos ( 



(12) 



and identical relations between [u' ,v') and (u,w). The 
idea now is to rewrite equations (jlip in terms of {x' , y') 
and (u', v'). As a way of example we have 



C Ux'x' 



+ S^U 



y'v' 



2csu 



x'y' , 



(13) 



where we have adopted the abbreviations sin 9 = s and 
cos 9 = c. Similar relations are obtained for the other 
second derivatives. Then u is replaced by m = cu' + sv' . 
These transformations allow one to rewrite the set of 
equations (fTTT) in terms of m^,-,, v'^,^, (with i',j' = x',y'), 
and ujj and v'^^. The transformed equations are listed 
in appendix [^ The two equations for uj^ and v'ff have 
now the form pH)) . not the shortest form (fTTj) . Invari- 
ance under rotation imposes that the form of the new set 
of equations should be identical to pT|) . and this leads 
to relations among the coefhcients. Since the invariance 
holds for any angle 9, it holds also in particular for spe- 
cial values of 9. For example, setting 9 — tt/2 simplifies 
greatly the algebraic equations (see appendix |X]) and this 
leads in a very simple way to the relations 



a' = 13, (3'=a, 7' - 7. 



(14) 



Then, imposing that the invariance holds for any angle 9 
leads to the additional relation 



7 = a-/3, 



(15) 



reducing the number of independent coefficients to two. 
Using these two relations allows to rewrite pT|) as 



pil = l3{Uxx + Uyy) -I- (a - l3){Uxx + Vxy) 



(16) 



PV = l3{Vxx + Vyy) -I- (a - l3){Vyy + Uxy), 

or in the more compact, vectorial form 

pu = /3V2u+(a-/3)V(V-u), (17) 

which is nothing but the classical elastodynamics equa- 
tions for an isotropic material. This concludes our proof 
for an isotropic 2D material. 

Equation ([T7| can be rewritten as 



pxii = djU. 



J "I J 



with 



2/3e,j- + (a - 2(i)ekk5^, 



(18) 



(19) 



which is the Hooke law (recall that e.^ = ^[di\ij + djVn) 
is the strain tensor). 

In conclusion, by evoking only elementary space sym- 
metries (and no notion of torque, or that the forces are 
of short range), we arrive to relation (|18|) . which implies 
the Cauchy principle, i.e. that the force is localized to 
surface only, and Hooke law ([TS]). In addition, the stress 
tensor (J19p is automatically symmetric. This section cor- 
responds to the main message of our paper. The rest is 
dedicated to extension to three dimensions and to other 
symmetries in two and three dimensions. 



III. EXTENSION TO OTHER 2D SYMMETRIES 

Let us still consider the 2D case, but other symme- 
tries, namely the hexagonal and square symmetries. The 
same transformed equations (see appendix]^ upon a ro- 
tation by an angle 9 can be used. For hexagonal sym- 
metry, we have to impose that Eqs. (|11|) are invariant 
for 9 = 7r/3. Simple algebraic manipulations lead to the 
conclusion that the six coefficients are related exactly by 
relations (fT4)) and (|15p . Therefore we get the same con- 
clusion as in the isotropic case. Actually it is well knowr^ 
in elasticity theory that in 2D the hexagonal case is iden- 
tical to the isotropic one (this ceases to be valid in 3D). 

For the square symmetry, invariance under a rotation 
7r/2 (already used in the isotropic case) has led to the 
relations (HH), entailing that there are three independent 
coefficients. The equations for the displacement field read 



pil = aUxX + PUyy + JVxy, 

pi) = Pvxx + avyy + ^Uxy 



(20) 



It is a simple matter to see that these equations can be 
put in the form ([T8| with 



a^xx + {1 - PYvy, 

yy OtfZyy + (^7 P J^XX 



^XX 
O 

^xy — ^yx ^P^xy 



Here again we recover the same conclusion: the symme- 
try of the stress tensor is automatically satisfied^ 



IV. EXTENSION TO THREE DIMENSIONS 

The extension of the above approach to three dimen- 
sions is straightforward and follows exactly the same 
spirit. Therefore we will keep this section brief. Let 
u,v,w denote the three components of the displacement 
field. Like in 2D there are no first derivatives entering the 
displacement equations. There are now 18 different sec- 
ond derivatives oiu,v, w. However several coefficients are 
zero. Indeed, using the 3D analogue of transformation (|S]) 
(which means imposing the additional invariance with 
respect to the symmetry operation z — >■ —z, w — > —w) 



easily allows to deduce the 3D analogue of (|TT|) 



pu 



fiUy 



IVxy 



5Uz 



PV = a'Vxx + fi'Vyy + "f'Uxy + 5'Vz 



(21) 



pw 



= a"'Wxx + l3"Wzz + l"uxz + &"Wyy + '-"v-. 



yyz 



There are a priori 15 different coefficients. For an 
isotropic material we can use the invariance with respect 
to a rotation with an arbitrary angle 9 and around an ar- 
bitrary axis (there are three independent axes). We first 
perform a rotation around the z-axis. The transforma- 
tion of (x, y) and (u, v) are identical to (J12p while z and 
w are unchanged. Transforming the (x, y, z) and (u, v, w) 
into (x' ,y',z') and [u' ,v' ,w') and writing the equations 
for it', i;' and ii)' and requiring invariance of ()21|) (see 
appendix |B| one straightforwardly obtains the following 
five relations 



/?' 



a, a 



/3, 7 = y, 5^5', 



(22) 



Using rotations around the two other axes {x and y) one 
obtains eight additional independent relations. In total 
we have thus 13 relations and we are left with only two 
independent coefficients in equations (I2ip . The reader 
can read off the full relations between coefficients just by 
confronting (PT|) with the following set 

pu = aUxx + PiUyy + Uzz) + {O- ~ P){Vxy + Wxz) 

pV ^ aVyy + fi{Vyy + Vzz) + (^ - (i){Uxy + W x z) (23) 

pw = aWzz + PiWyy + Wzz) + (« - fi){UxZ + Wxz) 

It is a simple matter (following exactly the same proce- 
dure as in 2D) to show that the vectorial form of this set 
is identical to (IT7)) and formally the relations (IT51) and 
P^ still hold (with now i,k ~ x,y,z). This concludes 
our analysis for a 3D isotropic material. 

It is easily seen in appendix [B] that for a cubic symme- 
try we are left with three independent coefficients, and 
the set of equations is similar to (|23p with the last coeffi- 
cient (a — P) replaced by 7. The equations can be put in 
the form (fTS]) , with the stress tensor having the following 
expressions for diagonal elements. 



(Jxx 


= atxx + (7 - 


- P)[iyy^tzz), 


'^VV 


= aCyy + (7 - 


- (3){exx + ezz), 


Ozz 


= aezz + (7 - 


- I3){exx + <^yy) 



and for non diagonal elements 



2/?e, 



(24) 



That is to say, the 3D cubic system equations can be 
derived fully from basic symmetry considerations. 



V. DISCUSSION AND CONCLUSION 

The main message of this article is that for high sym- 
metry crystals (isotropic, square and hexagonal symme- 
try in 2D, isotropic and cubic symmetry in 3D) the elas- 
ticity equations can be fully derived from just evoking 



basic symmetry considerations. From this approach, it 
follows automatically that (i) the force is a divergence of 
a second rank tensor (stating thus that the force acts on 
the surface only), and (ii) that the stress tensor is sym- 
metric. We find that this is interesting inasmuch as that 
these two points are usually motivated by physical prin- 
ciples: short range interaction for point (i) and absence 
of torque for point (ii). Therefore, the basic symmetries 
of these crystals induce naturally the symmetry of the 
stress tensor. 

Suppose we have a cubic crystal and consider a cu- 
bic element (in the continuum spirit). In textbooks, the 
force per unit surface acting on face "a;" (we mean the 
face having its normal along the a;-direction) and directed 
along y, is denoted as (Tyx, while that acting on face "y" 
and directed along x is denoted as <7xy Since we have 
exploited the basic spatial symmetries, the two faces are 
equivalent, and therefore the equality ayx = (Jxy, follows 
naturally and is encoded in the basic space symmetry. 
Therefore, evoking a zero torque is not only redundant, 
but also conveys the wrong impression that a physical 
principle would be needed. Had we considered a rect- 
angular symmetry (in the 2D case, for example), then 
no equivalence between the two faces is dictated by the 
spatial symmetry, and it seems necessary, at the present 
state of our understanding, that the zero torque condi- 
tion would be required. 

Note that research on elasticity theory with the aim of 
providing a deeper understanding has given rise to several 
papers. The most relevant in the spirit of the foundation 
of elasticity theory are the following two research. The 
work by Dimitryev who presented a nice study consisting 
in mapping the extended electrodynamics equations onto 
the theory of elasticityiS. This work has discussed the 
connection between the two theories, and entirely differs 
from our which raises the question of how to present the 
elasticity theory itself from elementary symmetries. The 
other noteworthy study is due to Borg^^ who presented a 
combination of dimensional analysis and symmetries and 
applied it to the derivation of the strain energy. How- 
ever, (i) the symmetry of the stress tensor is taken as a 
postulate, while here we have shown that it follows auto- 
matically from symmetries, (ii) he used the Hookes law, 
while here we have shown that it follows naturally from 
elementary space symmetries. 

Finally, we would like to point out a certain redun- 
dancy in some textbooks. Indeed, as soon as one assumes 
that the elastic energy is a function of the deformation 
tensor only (actually of its invariants), then this implies 
automatically point (i). Indeed, let the total elastic en- 
ergy reads 



fi^v)dV 



(25) 



where / is the energy density. A variation of the energy 
reads 

5F = J ^Se,,dV = i y ^{d,5u,+d,6ui)dV (26) 



where we have used the definition of the strain tensor eij 
in terms of the displacement vector. Since i, j are dummy 
indices, one can rewrite the last term as 



5F= f ^{d,duj)dV 



= - d,. 



where we have used the divergence theorem, and where 
Tij is the normal to the elastic surface body. Actually, 
if we have a finite body there are also surface forces gi 
(like the pressure field, for example) that produce a work, 
given by the force times the displacement that has to be 
added above, — J gjSujdA. 

The first term on the right hand side of (l?f)) is a prod- 
uct of some entity, di{jr~), times a displacement, but 
this is nothing but the very definition of the work done 
by internal stresses. In other words i9i(g^) is the force, 
which naturally reads as the divergence of a second or- 
der tensor aij = -fp—- At equilibrium the total variation 

vanishes. Thus, the bulk term in p7|) vanishes (corre- 
sponding to V<T = 0) as well as the surface term, leading 
to gi = (JijUj, which is nothing but the boundary condi- 
tion at the external surface. 

Acknowledgement: CM. Acknowledges financial sup- 
port from CNES (Centre National d'Etudes Spatiales) 
and ESA (European Space Agency). 



Appendix A: The 2D case 

Using the transformation ([T^ the set of equations PTj) 
transforms into 

p(cu + sv) =Au^,^, + Bu , , -t- Cv'^,^, 



+ ^<'y' + ^<'y' + ^<'v' 



p{cv' - su') =A'u'^,,^, + B'u'^, + C'v'^,^, 



(Al) 



where 



y y 



A ^ (ac^ +/3s2)c + 7cs^, 
B ^ (as'^ +(3c^)c-js'^c 
C = {ac^ +/3s^)s-jc'^s, 
D = (as^ + /3c^)s -I- 7c^s 

F ^2{a- P)cs^ + c{c^ - 5^)7 
A' = -{a'c^ +P's^)s-ic^s, 
B' = -{a's^+/3'c^)s + 'y'c^s, 
C' = (a'c2+/3's2)c-7'cs2, 
D' = {a's^+(3'c^)c + 'y'cs^, 
E' = 2(/3' - a')c.s^ + c{^ - 5^)7', 
F' = 2{a' - I3')c's + sic^ - s'W 



Multiplying the first (resp. second) equation in (jAip by 
c = coaO (resp. s = sin0) and the second by s (resp. c) 
and making the difference (resp. the sum), one obtains 
the equations for il' and v': 



y'y' 



pu' = {cA - sA')u'^,^, + {cB - sB' 

+ {cC - sC')v',,,, + {cD - sD')v'y,y, 
+ [cE ~ sE')u'^,y, + {cF - sF')v',,y, 

pv' = {sA + cA')u'^,^, + {sB + cB')u'y,y, 
+ {sC + cC')v',,,, + {sD + cD')v'y,y, 
+ [sE -f cE')u'^,y, + {sF + cF'W, 



(A2) 



,yl 

Comparison with (jlip tells us that the invariance is pre- 
served if and only and if we have the following relations 

cA - sA' = a, cB- sB' = (3, cF ~ sF' = 7, 

cC - sC' = 0, cD- sD', cE- sE' ^0 

sC + cC' = a', sD + cD':^l3', sE + cE' ^ j' , ^ ' 

sA + cA' = 0, sB + cB' = 0, sF + cF' = 

For a square (as well as an isotropic) crystal this relations 
holds for 9 — tt/2. The algebra is quite simple (since c — 
0, s = 1), and we immediately obtain (|14p. This means 
that for a square symmetry only relation (fT4|) holds. For 
an isotropic material we still can impose that the set 
(IA3|) holds for any 9 (by using the already obtained result 
(TTil) . which greatly simplifies the algebra). This yields 
relation (|15p . Finally for an hexagonal symmetry, we 
must impose invariance with respect to 9 = vr/S (c = 1/2 
and s = •\/3/2). Reporting into (jA3| leads to quite simple 
algebraic equations leading to relations (fni) and (|15p . 
In other words, the hexagonal symmetry obeys (in two 
dimensions) exactly the same set of equations as isotropic 
materials, a well known fact^. 



Appendix B: The 3D case 

The 3D case presents no special challenge, and the rea- 
soning is exactly the same as in 2D. We start from (|2T|) 
and performs a rotation around the z-axis. In this case 
(x', y\ z') = {ex — sy, sx + cy, z) and identical relations 
between the displacement components. The transformed 
equation in the new coordinate system reads 



p(cu' 



= r, 



2D 



- <>(.CU'^'Z' 



') 



+ l{cw'^j^, +sw' ,^,) 



p{cv' -su') ^T^jj + S'{c 



pw 



+ l-'iCWy,,, - SW'^,,,) 

+ p"w'^,^, -f i'{c^u'^,,, + csv'^,,, + csu'y,^, + s^v'y,^,) 

+ ''"{c^Vy,^, - CSV'^,^, - CSu'y,^, + S^u'^,^,) 

+ '5"(S2<,,, + C'w'y,y, - 2CSW',,,y,) 

(Bl) 



where T2D and Tl^jj are the right hand sides of the first 
and second equations in (|A1I) . respectively. The next 
step is to perform the same manipulations done for the 
2D case (i.e. extract equations for u' and v'; that for w' 
is already explicit in the set (IBljl ). Imposing invariance 
of obtained equations (i.e. that they must have the form 
(PTjl ) yields the same relations ((Ti| and ([T5|) (obtained in 
2D) plus two new independent relations (implied by the 
third variable w): 



6^6', i = i'. 



(B2) 



Performing rotations about the other two axes x and y 
axes yield new additional relations among coefficients. 
Actually, it is even much simpler to deduce the new rela- 
tions, just by successively exchanging the role of w and u 
and then the role of w and v (without making any alge- 



braic calculation) . For example the above equality 5 — 5' 
means (when inspecting the set (PTjl ) that the coefficient 
Uzz in the u-equation is equal to the coefficient of Vzz 
in the u-equation. This obviously implies that the coef- 
ficient of Wxx in the w-equation should be equal to the 
coefficient of Vxx of the w-equation. For the same rea- 
son, the coefficient of Wyy in the w-equation should be 
equal to the coefficient of Uyy of the u-equation. Simi- 
lar reasoning can be made for the relation l = i' . These 
considerations lead to 



5" = 5 = /3, 



i = i = 7, 7 =7. 



(B3) 



Using the relations dUl), ^^, (lB2| and ([53]), leads to the 
final system p3|) . For the cubic case, the relation p^ 
(which follows by imposing an arbitrary rotation angle) 
does not hold and the equations are obtained by replacing 
the last coefficients {a — f3) by 7. 
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L. D. Landau and E. M. Lifshitz, Theory of elasticity (El- 
sevier, 2007). 

^ R. P. Feynman, R. B. Leighton, M. Sands, Feynman lec- 
tures on Physics: Mainly electromagnetism and matter 
(Addison- Wesley, 1964). 

^ A. E. H. Love, A treatise on the mathematical theory 
of elasticity (Dover Publications, 1944). In p. 76-78 the 
changes of momentum is defined, and the torque is as- 
sumed to vanish leading to the symmetry of the stress ten- 



* S. Timoshenko and J. N. Goodier, Theory of elasticity (2nd 
edition, McGraw-Hill, 1951) At p. 4, authors find the sym- 
metry of the stress tensor from the equation of equilibrium 
(no torque). 

^ C. Kittel, Introduction to solid state physics (7th edition, 
John Wiley, 1996). At p. 82, author writes: The number 



of independent stress components is reduced from nine to 

six by applying to an elementary cube the condition that 

the angular acceleration vanish, and hence that the total 

torque must be zero. 
® For ease of notations, Ui will indicate a spatial derivative, 

see the next Section. 
^ The more common Young modulus E and Poisson ratio 

v are algebraic combinations of A and /j.: E = /i(3A + 

2/x)/(A + n) and u = A/2(A + fi). 
* We are assuming here that forces are local. 
® More generally, we can show that Eqs. pip can always 

be written in the form (|18|) . However, the resulting stress 

tensor is symmetric if and only if a' — p. 
^^ V.P. Dimitryev "Electrodynamics and Elasticity" Am. J. 

Phys.71 (9) 952-953 (2003). 
^^ S.F. Borg "On an Application of Dimensional Analysis" 

Am. J. Phys.19 (2) 69-73 (1951). 



